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The nature of multiple electronic orders in skutterudite PrFe4Pi2 is discussed on the basis 
of a model with antiferro-quadrupole (AFQ) interaction of T3 symmetry. The high-field phase 
can be reproduced qualitatively provided (i) ferro-type interactions are introduced between 
the dipoles as well as between the octupoles of localized /-electrons, and (ii) separation is 
vanishingly small between the ri-T 4 1 ' ) crystalline electric field (CEF) levels. The high-field 
phase can have either the same ordering vector q = (1,0,0) as in the low-field phase, or a 
different one q = depending on the parameters. In the latter case, distortion of the crystal 
perpendicular to the (111) axis is predicted. The corresponding anomaly in elastic constants 
should also appear. The electrical resistivity is calculated with account of scattering within 
the CEF quasi-quartet. It is found that the resistivity as a function of the direction of 
magnetic field shows a sharp maximum around the (111) axis at low temperatures because 
of the level crossing. 
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1. Introduction 

Rare earth filled skutterudites are intensively studied in the last few years because of 
their variegated and complex behavior. In this paper we focus among others on PrFe4Pi2 
which exhibits a phase transition at 6.5K. A fascinating new experimental result on PrFe4Pi2 
is the discovery of a high- field ordered phase sharply located around the (111) direction of 
the magnetic field. 1 It is natural to expect that this new phase is related to the strange 
angle dependence of the electrical resistivity which shows a sharp peak around the (111) field 
direction. 2 The purpose of this paper is to provide a consistent account of the phase diagram 
of PrFe4Pi2 including both low- and high-field orders. On the basis of theoretical analysis, we 
predict possible lattice distortion, which leads to a symmetry lower than trigonal along (111), 
and elastic anomaly on entering the high field phase. 

Pr ions in the PrFe4Pi2 filled skutterudite have 4/ 2 configuration with total angular 
momentum J = 4. Under the point-group symmetry T^, the ninefold degenerate J = 4 
level is split into CEF levels. The nature of the low energy CEF levels is essential to the 
physical properties and electronic orderings of the system. The CEF splittings in related 
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systems PrRu4Pi2 and PrOs4Sbi2 are clearly known by neutron diffraction experiments. 3 ' 4 
However, the sequence of CEF levels in PrFe4Pi2 has not yet been established experimentally. 
Interpretation of neutron scattering spectra of PrFe4Pi2 remains ambiguous: polycrystalline 
samples do not show clear peaks and the peaks appearing in the single crystalline samples 
become clear only below the transition temperature. 5 The first attempt on determining the 
CEF scheme was the analysis of the magnetization anisotropy in the disordered phase. 6 It was 
found that the anisotropy can be reproduced either by schemes with ri— T4, T1—T5 or T3-T4 
low- lying levels within the point group Oh- Measurements of elastic constant found softening 
of the modes C\\ and C\\ — C12 towards low temperatures. The results were interpreted in 
terms of the T3 doublet ground state. 

We find that the T±— scheme in magnetic fields gives rise to a level crossing only for 
field direction (111). Our result generalizes the observation by Tayama et al., 1 who worked 
with the model with the point group Oh- We thus consider that the relevant CEF scheme in 
PrFe4Pi2 is the Y\-Y^ singlet-triplet scheme, which is in line with ref.8. This level crossing is 
the key aspect of the high-field phase and will be fully discussed in Section 3. If one takes the 
doublet ground state, on the other hand, the energy levels in magnetic field are rather different. 
With the doublet ground state, it is difficult to explain the high-field phase in terms of the 
level crossing. Analyzing the quadrupole pattern in the ordered phase at low temperatures, 
we shall clarify that a two-sublattice AFQ order of the T3 quadrupoles within the Ti-r^ 
subspace has a unique feature that does not share with a quadrupolar order in other CEF 
states such as the T3 doublet or the T$ quartet. 

This paper is organized as follows. In Section 2 single-site properties of our model is 
discussed. In Section 3 we examine the T3 quadrupole interaction model for the field direction 
(111) in detail. The effects of ferro-type dipolar and octupolar interactions on the phase 
boundary are analyzed. We also derive the Curie- Weiss temperature. In Section 4 the electrical 
resistivity is calculated using the quasi-quartet crystal field states of T function of 

H||(lll), and the angle dependence around the (111) direction is derived. In Section 5 some 
basic features of the T3 quadrupolar ordered phase in the triplet is discussed. Section 6 is the 
summary and discussion of our results. 

2. Singlet-Triplet CEF states and their splitting in magnetic field 

2. 1 Crystal-field states and T% quadrupoles 
We start with the single-site Hamiltonian: 

H ss = Hcef ~ <?/U B H • J, (2.1) 

where the second term represents the Zeeman energy with the angular momentum J, and 
the first term is the CEF potential under the tetrahedral crystal field. The CEF term can be 
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Hcef = W (x0 4 + (1 - x)O c 6 + yO\) , (2.2) 

where we have used the normalization convention 0\ = (j\ + Ji)/120 with J± = J x ± iJ y in 
4 = 0° + 50|. 

The term yOg has no effect on the cubic Ti and T3 levels, but mixes the cubic T4 and 
T5 triplet states. Then the new triplets T^ 1 -* and arise. The triplet can be written in 
terms of Oh triplets r^ "* and as 



ri 1} ,m) = d\rf\m) + \/i-d 2 |rf ,m) , 



,(0) 



(0) 
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We further introduce the notations for future convenience 

|r«o) = |o), |rW,±) = |±). 

The energies of the states of the pseudo-quartet are given by 



(2.6) 



(2.7) 



2W 



x - 4 + 2V(3 + 2x) 2 + 1008y 2 



S(ri) = W [108a; - 80] . 



(2.8) 



We note that the parameter y is indispensable for getting the T\ and states closely located 
to each other without the proximity of the T3 doublet state. 

Now, let us consider the 1^3 quadrupole moments under the tetrahedral symmetry T^. As 
an alternative to the pseudo-spin representation, 10 we directly compare matrix elements of 
the quadrupolar operators O = [3 J 2 — J (J + l)]/\/3 and 0\ = J 2 — Jy with respect to 
the tetrahedral triplet states {|ai), \cx2), |<S3)} = {|0), |+), |— )}, and the cubic T4 triplet 
states \di) (i = 1, 2, 3). Then we define the 3x3 matrices O2 and C% with fx = 0, 2 by 



(0%)* = {Oi\0$\a k ), (C%) ik = {(H\0$\a k ). 



(2.9) 
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C>2 = ciqC® + , 

O 2 = a Q C 2 -b Q C° 2 , (2.10) 

where a Q = (9d 2 - 2)/7 and b Q = -d^3{l - d 2 )/7. The limits d = 1 and mean the cubic 
triplets T4 and T^, respectively. Namely we recover 0\ = C® and 0\ = C 2 from (2.10) because 
of aQ = 1 and 6g = 0. The inverse transformation of equations (2.10) gives 



° 2 ~ J^TbJ)° 2 (4T65)° 2 



C 2 = 7At02 + 7At05. (2.H) 
(a 2 +b 2 ) {a 2 +b 2 " ' 



We remark that 0% with fi = 0, 2 can be diagonalized simultaneously by the basis set |±), |0). 
The eigenvalues can be obtained from the eigenvalues — 14/\/3(l, 1, —2) of C\, together with 
eigenvalues ±14, of C 2 by proper combination. Using equations (2.11) we obtain 

{Ol) 2 + (O 2 ) 2 = (a 2 Q + b 2 Q Wl) 2 + (C 2 ) 2 } = 7 -f{a% + b 2 Q ), (2.12) 

which is a scalar matrix in consistency with the fact that (Of>) 2 + (0|) 2 ^ s a Casimir operator. 
The factor a 2 Q + b 2 Q = (4 - Ihd 2 + 60d 4 )/49 reduces to unity with d = 1, and to 4/49 with 
d = 0. Namely the r^ -* triplet has a smaller quadrupole moment as compared with the 
triplet. 

2.2 Constraints on the singlet-triplet CEF splitting 

Analyzing the energy spectra of the Hamiltonian Tt ss = TtcEF — ^bH • J for different 
crystal field parameters, we can easily find that only the Fi-T^ 1 - 1 low-lying scheme gives a 
crossing point uniquely for the (111) direction of the magnetic field. The mixing parameter d 
has to be larger than d cr = \/l5/6 w 0.645, because for d < d cr a level crossing occurs also 
for H||(001). This conclusion comes from the comparison of the inter- and intra-level matrix 
elements of J z within the Fi-r^ quasi-quartet. These are ±2\/l5<i/3 and ±(5/2 — 2d 2 ), 
respectively. 

Furthermore, the existence of the high-field phase imposes a strong constraint on the 
singlet-triplet CEF splitting. With the singlet CEF ground state, the quadrupolar order must 
be an interaction induced order. We need a finite coupling constant between the quadrupoles 
for the phase transition even at zero temperature. In the case of a large energy separation A, 
on the other hand, the transition changes to first-order from second-order for A > A ter , where 
A tcr means the tricritical value of the energy separation. Loosely speaking, in a mean-field 
theory we expect that A tcr has the same order of magnitude as the transition temperature 
Tq = 6.5K. Experimentally the phase transition at Tq is second-order, hence we obtain the 
upper bound of the CEF splitting of about 7K within our model. On the other hand, the high- 
field phase occurs with infinitesimal interaction because of the level crossing. If the interaction 
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parameters responsible for the low field AFQ order works also in the high-field phase, the 
observed phase diagram suggests very small splitting of the quasi-quartet. Otherwise, the 
large interaction necessary for the induced AFQ order makes the high-field phase much wider 
than actually is. The critical field H cr becomes also too large. 

Therefore in the following we assume a singlet-triplet energy gap vanishingly small. This 
assumption E(T^) — E(T±) = defines a relation between the tetrahedral crystal field pa- 
rameters x and y as given by eq.(2.8). Together with the constraint on d discussed earlier, 
we choose d = 0.9 which gives x = 0.98 and y = —0.197, and we use W = — 0.9K. Figure 1 
shows evolution of the crystal field levels as the magnetic field increases with these parameters. 
Increasing the strength of the magnetic field, a level crossing occurs in the ground state. We 




note that the behavior remains similar in the presence of small positive A. 

We also derived the field evolution of the CEF states for field directions (001) and (110), 
which are shown in Fig. 2. We can observe that no level crossing occurs in the ground state for 
these field directions. Therefore, the effect of application the magnetic field in these directions 
is expected to be the suppression of the antiferro-quadrupolar ordered phase. 
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3. Phase diagram in the multipolar interaction model 

3. 1 Instability of the paramagnetic phase 

We shall work with the minimal model for intersite interactions which can reproduce the 
two ordered phases at least qualitatively. We include such octupoles that reduce to T^, u in 



the cubic symmetry. These are given by T§ = (J x Jy — J|J x )/6, Ty = (JyJ? — J x J y )/6 and 
Tz = (J z Jx ~ JyJz)/®, where the bar means symmetrization. Then we consider the following 
model with the nearest-neighbor intersite interaction: 

Hint = \Y. [ X ^°h ■ O 2:j + Ol t • 2 2J ) + A oct Tf • Tf + X dtp 3, t • J,] , (3.1) 

with Ar 3 > and X oc t < 0, Adi p < 0. Namely we take antiferro-type quadrupolar and ferro-type 
dipolar and octupolar interactions. The dipole and 1^3 quadrupolar interactions are consistent 
with the invariant form allowed by symmetry of the bcc lattice. 11 On the other hand, for 
simplicity we have neglected the mixing term J • T' 3 = J X T^ + J y Ty + J z Tz between the 
dipoles and T@ octupoles, although the mixing is allowed by the tetrahedral symmetry. In 
the following we constrain the magnitudes of parameters such as the coupling constants and 
the CEF splitting. 

We take the model given by 7i ss + Timt , and obtain the phase boundary in the plane of 
magnetic field and temperature. The nine crystal field levels labeled by \l) are derived for 
different values of magnetic field by solving H ss \l) = Ei\l). In deriving the phase boundary, we 
keep only the lowest six levels for simplicity of the numerical calculation. As seen from Fig. 1, 
the three higher levels neglected here repel with lower levels, and go up in energy. Thus their 
neglect will not influence the low-temperature behavior significantly. 
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From the six levels labelled by \l) with I = 1, 6, we have 35 (= 6 x 6 — 1) independent 
pairs k, I describing multipole operators at each site. We have 15 symmetric combinations 
of the pairs, and arrange them as increasing order of lOfe + / with k < I. Then we define 
X a = \k)(l\ + \l){k\ w hh a = 1, 15 according to the ordering of the pair. Similarly we define 
15 antisymmetric combinations X@ = i • (|A;)(Z| — with (5 = 16, ...,30. The remaining 

five operators are diagonal in k and I. The general form of the multipolar interaction is given 
by 

= ~ \ E E V f X ? X i = ~ \ E ^ ■ % ■ X, , (3.2) 

where X = [X 1 , X 2 , .., X 35 ], and i,j are site indices. The susceptibility matrix can be ex- 
pressed as 

X = (l-Xo-v)' 1 -Xo, (3.3) 
where xo is the single-site susceptibility matrix whose elements are given by 

( X0 r = E p{Ek ^z p E {El \ k\x«m\xP\k) , (3.4) 
k,i 1 h 

where p(Ek) = exp(— fiE^)/ exp(— f3Ek). In the right-hand side of eq.(3.4), we only need 
terms with k ^ I, since the interaction Hamiltonian (3.1) does not have the diagonal part of 
operators X a . The calculation of (xo) a ^ c & n easily be performed with use of the magnetic 
eigenstates, which are nondegenerate except at the level crossing point. It is important here 
to include renormalization of the external magnetic field h = qixbH. Namely we consider the 
effective magnetic field given by h e s = h — zXdip{Jin) where (Jm) depends on the effective 
field h c g and temperature T. The effective field should be determined self-consistently. 12 
Actually we first take h e s as a given field and derive the matrix xo under h e s. To each value 
of h e ff, we determine h by the relation h = h e s + z\di p (Jm)- 

After Fourier transformation we obtain the susceptibility matrix as a function of q: 

x(q) = (l-Xo-^(q)) _1 -Xo, (3.5) 

where V al3 (q) is the Fourier transform of V^f ■ The interaction matrix V a/3 (q) can be classified 
in a form 

F^ (q) = vg (q) + Cf (q) + <(q)- (3-6) 
A second-order phase transition is characterized by the divergence of the susceptibility matrix. 
This is equivalent with condition 

det (l - xo • nq)) =0, (3.7) 

where 1 means the 35-dimensional identity matrix. The solution of eq. (3.7) gives the instability 
condition of the paramagnetic phase. Note that the phase boundary may not coincide with 
the instability line if there is a first-order transition. In this paper, we do not go into detailed 
inspection about the order of the transitions. 
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3.2 Ordering vector and instability lines in magnetic field 

Under the reduced symmetry with finite magnetic field along the (111) direction, the 
original Th system allows only two different order parameter symmetries: the C^iTi) singlet 
and the Cs(T^) doublet representations. 10 Within the two-dimensional local Hilbert space, 
which is formed by the closest crystal field states near the level crossing point, the order 
parameters with C^iT^) symmetry have nonzero inter-level matrix elements. If the high-field 
phase has the same ordering vector as the low-field AFQ phase, both phases can be connected 
smoothly. Alternatively, the high-field phase can have a different ordering vector q without 
the smooth connection. 

We derive the instability toward the low-field and the high-field phases from high tem- 
peratures with several sets of interaction parameters. We study both cases with the ordering 
vector q = (1,0,0), which leads to a two-sublattice order on the bcc. lattice, and q = 0. 
The quadrupolar coupling constant Ar 3 is determined by the zero field transition temperature 
Tq = 6.5K. In order to determine this quadrupolar coupling constant within the lowest six 
levels, let us consider the mean-field form of the Hamiltonian H ss + Hint at H = 0: 

Hmf(Ti - ri 1} - r 3 ) = a l r 3><r?l + z\ r3 ((o° 2 ) A{B) o 2 + (o 2 2 ) A{B) o 2 2 ) , (3.8) 

n=l,2 

where the quasi-quartet level is taken as the origin of energy, and the energy of the 1^ level is 
A. We perform the Landau expansion of the free energy. The quadrupolar phase transition, 
which is second order at H = 0, is determined by zero coefficient a q of the quadratic term. 
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Fig. 4. Temperature-magnetic field phase diagram with Ar 3 = 8.93mK, Xdip — — 102. 5mK, and 
A oc t = — 11.39mK. Dashed line shows the observed phase boundary. 



Two different solutions can be obtained depending on the magnitude of ferro-type coupling 
constants: (a) both low- and high-field phases have the same ordering vector q = (1, 0, 0); (b) 
the high-field phase has the ordering vector q = 0. In the case (a) we obtain either continuous 
phase boundary or two separated phases (see Fig. 3). For example, with X^p = — 125mK 
and X oc t = — 8.55mK we obtain the continuous behavior. While with the same Xdip but 
Xoct = — 8.625mK, two phases are separated. An example for case (b) can be seen in Fig. 4, 
where the parameter are A^ p = — 102. 5mK, X oct = — 11.39mK. 

The effects of the ferro-type interactions on the phase boundary depend on the multipolar 
matrix elements and the quadrupolar-dipolar-octupolar mixing due to the symmetry lowering 
as we increase the magnetic field. We found that the dipolar matrix elements are effective 
in the reduction of the transition temperature in the low- field regime (H < H cross ), while 
the T@ octupolar matrix elements have larger effects in the high-field regime (H ~ H cross ). 
This latter behavior arises because the inter-level octupolar matrix elements between T^-I^ 
levels become significant as we increase the magnetic field. In the case of larger ferro-octupolar 
coupling, the high-field phase thus takes q = as the ordering vector. 
3.3 Magnetic susceptibility and Curie-Weiss temperature 

We examine whether the parameters used in the previous calculation are consistent with 
other observables. For this purpose we discuss the temperature dependence of the magnetic 
susceptibility in the high-temperature phase. Comparsion of the Curie- Weiss temperature 



We obtain 



1 2 4 [exp(-A/r)(16A - MOT) + U7(a 2 Q + b 2 Q )A + MOT] 

■q- 2 *Ar 3 z A r33 TA(4 + 2exp(-A/T)) 



(3.9) 



With the parameters used in Fig. 1 we obtain Ap 3 = 8.93mK from a q = 0. 
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between theory and experiment will provide a test of our model. The positive Weiss temper- 
ature found experimentally is the most straightforward evidence of the ferro-type interaction 
between the dipoles. We remark that a ferromagnetic phase transition has been found in 
Pr :E Lai_ a ;Fe4Pi2 for x = 0.85, 13 which also indicates the strong ferro-dipolar interaction. We 
have further introduced ferro-type interaction between the T@ octupoles. In the tetrahedral 
symmetry, these octupoles have the same symmetry as the dipoles (r4 U ), therefore they cannot 
be distinguished from the dipoles experimentally. 

We start with the single-site susceptibility matrix. Instead of X a in eq.(3.4) as multipole 
operators, we can equivalently use another set of operators with Z 1 = J z and Z 2 = Tz . 
We obtain the 2x2 matrix for the single-site susceptibility as 



(x r = £p(£* 



(3.10) 

k 

The susceptibility in eq.(3.10) has the Curie term in contrast with eq.(3.4). Then the magnetic 
susceptibility xm including the intersite interactions is given by 



Xm : 

where the matrix V has the form 



(3.11) 

J li 




I | "'" ° ] , (3.12) 

zXoct J 

with z = 8. Because of the off-diagonal component Xq 2 an d Xo 1 ' octupolar interaction influ- 
ences the magnetic susceptibility. The susceptibility given by (3.11) diverges at the Curie- Weiss 
temperature. This condition is equivalent with the divergence of the q = component of the 
susceptibility (3.5) in the absence of external magnetic field. 

With the set of parameters used in the left panel of Fig. 3, namely A^p = — 125mK and 
^oct = — 8.55mK, we obtain 6* = 4.41K. On the other hand, another set A^p = — 102. 5mK and 
Kct = — 11.39mK, which is used in Fig. 4, gives 6* = 4.12K for the Curie- Weiss temperature. 
These are to be compared with the experimental value 0* ~ 3.6K. 6 

4. Electrical resistivity in magnetic field 

In this Section we analyze the effect of the level crossing on the electrical resistivity. 
We are especially interested in the sharp enhancement of resistivity around the (111) field 
direction. We use a very simple model where conduction electrons are scattered by magnetic 
exchange and aspherical Coulomb interaction. The change of population of the crystal field 
levels causes a temperature and magnetic field dependence in the resistivity. We follow the 
method presented in refs.14, 15, and calculate the resistivity in the disordered phase. 

Using the notation of ref.16, the contribution to the resistivity from the CEF effects can 
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Fig. 5. Left: Electrical resistivity calculated as a function of H|| (111) at T = 2.3K and 3K using the 
crystal field scheme shown in Fig. 1. Right: Electrical resistivity calculated as a function of angle 
9 (sec text) near the (111) direction at H = 25. 3T for T = 2.2K and 3K. In both figures the 
resistivity is in arbitrary unit and r = 1/2. 



be described as 

PCEF = Po[rTr(PQ M ) + (1 - r)Tr(PQ A )] , (4.1) 

where r means the ratio between the two scattering terms. The matrix P has the form 

_ P (E k )f3(E k - EQ _ ^pjEMEjMEk - E t ) 
m l _ e -p{E k -E l ) 2cosh\(3{E k -E l ) K ' 

in the notation of eq.(3.4). The matrices Q M and Q A represent the magnetic exchange and 
aspherical Coulomb scattering, respectively. The latter comes from the quadrupolar charge 
distribution of the Pr 3+ ions. They are given by 

Oh = \(k\JM 2 + \(k\Jy\l)\ 2 + \(k\J z \l)\ 2 , 

Qti = E \( k \ Y 2 m \l)\ 2 , (4-3) 

m=-2 

where Y™ are the Stevens operator equavilents of the spherical harmonics for L = 2. 

With the crystal field parameters and the corresponding level scheme presented in Fig. 1, 
we calculate the electrical resistivity considering the lowest six levels as a function of H|| (111). 
We used the dipole interaction parameter A^ p = — 102. 5mK of Fig. 4 for the determination of 
the external field (H) from the molecular field. The octupolar interaction does not contribute 
to the molecular field in our model of eq.(3.1). The result is shown in the left part of Fig. 5. 
Due to the level crossing at H ~ 23T, the resistivity shows a corresponding peak at low tem- 
peratures. This peak can be interpreted as an enhanced contribution to P k i from degenerate 
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states. 14 Let \k) and |/) be the ground and the first excited states, respectively. If higher levels 
can be neglected, we obtain 

P(E k - E t ) 



This function shows a peak around E k — E\ = 0, which means an increase in the resistivity. 

We calculated also the electrical resistivity as a function of the polar angle 9 of the mag- 
netic field. This 9 describes the rotation of the field from H||(001) to H||(110) within the 
[1,-1,0] plane. The right panel of Fig. 5 shows the 9 dependence of the calculated electrical 
resistivity for a field close to the level crossing point. The resistivity is enhanced around the 
(111) direction (9 ~ 55°), and the peak becomes sharper as temperature decreases. In this cal- 
culation we have taken r = 1/2. We have checked that the results for H- and ^-dependences 
hardly change if we take different values of r. Experimentally, an even sharper resistivity 
peak has been found at T = 0.36K and H = 8T. 2 However, according to our model, we are 
already inside the ordered phase at this T and H as shown in Fig. 4. Hence our approxima- 
tion assuming independent impurities cannot be applied for quantitative comparison with the 
experiment. 

5. Quadrupole patterns at low temperatures 

In Section 3 we derived the instability of the paramagnetic phase without specifying the 
AFQ pattern, since the transition temperature is independent of the pattern as long as they 
belong to the T% order-parameter space. However, at low-temperatures the stability of the 
ordered phase should depend on the AFQ pattern. Hence it is interesting to study the basic 
features of the T% quadrupoles in the ordered phase. For simplicity, we examine the problem 
in the ground state and take H = 0. 

First, let us consider a two-site problem with sites A and B with the interaction Hamil- 
tonian within the subspace: 

W(r 3 ) = Ar 3 (O 2 %0 2 % + 0\ a O\ b \ (5.1) 

We take the nine-fold basis \k)A\l)B for the two sites, where k,l can be 0, + or — . By diag- 
onalization of H(Ts) we find that the nine-fold degeneracy splits into a six- and a three-fold 
multiplets. With antiferro-type interaction (Ar 3 > 0), the ground state is the six-fold multiplet 
with states 



-)a\0)b, 

-)a\+)b), |6) = ^(|+)a|+>b-|->a|->b). (5.2) 
To obtain insight into the nature of degenerate wave functions, we associate the states of the 
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= \0)a\+)b, 
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|0> ~ Y 2 ° ~ 3z 2 - r 2 , |+) ~y 2 2 ~ (x + iy) 2 , \-) ~Y 2 2 ~ (x-iyf. (5.3) 

If we choose the state \0)a, namely (3z 2 — r 2 )^, the degeneracy of 1 1) , |2) allows combinations 
|1) ± 1 2) to make (3z 2 — r 2 )A{x 2 — v 2 )b and (3z 2 — r 2 ) A (xy)B as the degenerate ground 
states. Similarly, combinations |5) ± |6) give the other ground states as (xy)A{x 2 — y 2 )B and 
(x 2 — y 2 ) A {xy)B- We remark that the degeneracy 6 comes from the number of ways for choice 
of different orbitals at two sites. 

The states |0), (|+) + \—))/y/2 and (|+) — \—))/y/2 are simultaneous eigenstates of O 2 
and 2 . The degeneracy of the same nature can also be found in the mean-field theory. In 
fact, the two-site mean-field theory at zero temperature becomes exact because of the absence 
of off-diagonal elements in the Hamiltonian with respect to the basis given by eq.(5.2). This 
situation remains the same in the bcc lattice, which can be separated into A and B sublattices. 
Hence with only the I^ intersite interaction, the ground state of the bcc lattice is degenerate. 
Recognizing this situation, it is instructive to consider the mean field of H{Pz) at site B: 

H B = Ar 3 ({O 2 ) A O° 2 + (0 2 ) A 2 ) = Ar 3 (Q A O 2 + q A 2 ) , (5.4) 

where we have used the notations {0 2 )a = Qa and {0 2 ) A = q A - If we choose the state (0)^, 
TCb is diagonal with respect to the basis set {|q;i), \ot2), = {|0), |+), | — )}• Namely we 

obtain the matrix representing the mean-field at B site: 



hBQ = t — {{ak\H B \ai)} = 
Ar 3 





/ 


2 


\ 


2c 




-1 












- 1 / 



(5.5) 



where c = (l4/\/3) (oq + bq) and we have used Q A = 28/ VSciq and q A = — 28/\/36q. We 
note in eq.(5.5) that the splitting into a singlet and a doublet is due to the mixing of moments 
O 2 and 2 in the molecular field. This quadrupolar field can be expressed as 

QaO% + q A 0\ = 28/^3(4 + b 2 Q )C Q 2 ~ 3z 2 - r 2 . (5.6) 

The degeneracy in hso with respect to the states |+) and |— ) is also indicated in the two-site 
calculation as seen in eq.(5.2). 

Another choice for the states on A is (|+) ± |— }) A . With this choice we get the following 
mean field at B site: 

/ -2 



hB± 



V 








1 

±3 



\ 

±3 

1 / 



(5.7) 



where the basis is the same as before: ||0), |+), |— )}, and we have used Q A = —lA/^/3aQ — 
146q, q A = 14/v / 3&q — 14oq. After diagonalization of hs±, we have the same eigenvalues with 
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those in hso of eq.(5.5). This means a degeneracy with respect to the states \0)b and one of 
(|+) =p \ —))b- The quadrupolar fields coming from the sublattice A are now given by 

Q A O 2 + q A O 2 = -U/V3(a 2 Q + b 2 Q )(C° 2 ±V3C 2 )~3y 2 -r 2 , 3x 2 -r 2 , (5.8) 

respectively. Even in the cubic limit, pure staggered order of 0\ does not exist, and the 
homogenous component of 2 is always induced. This is interpreted by the Landau free energy 
expansion in terms of order parameter components that the third-order term (T3 ® Yj, ® Y-$) 
is present where one of them is homogeneous, and two of them are staggered. 

Thus, we conclude that a model with nearest-neighbor AFQ interaction of ]?3 quadrupoles 
has a degeneracy with respect to different kinds of quadrupolar patterns. This degeneracy 
comes from the situation that we cannot construct an AFQ order with "anti-parallel" ar- 
rangement of the quadrupolar moments within the triplet, contrary to the cases of classical 
s = 1/2 spins, and AFQ pseudo spins within the Y3 doublet state. In other words, we have to 
choose an orthogonal arrangement. Choosing 3z 2 — r 2 -type quadrupole on one site, which can 
be realized by the state |0), there remains two possibilities for choosing an orthogonal state on 
the other sublattice. Namely, 3y 2 — r 2 and 3x 2 — r 2 are equivalent choices. In this sense, the T3 
antiferro-quadrupolar ordering model within the triplet state is similar to the three-state anti- 
ferromagnetic (AFM) Potts model, 17 which possesses a macroscopic degeneracy in the ground 
state. The degeneracy in our case is broken by dipolar and/or octupolar intersite interactions. 
With applied magnetic field, on the other hand, the Zeeman energy breaks the degeneracy 
even without magnetic interactions. If a weak magnetic field is applied along the direction 
(001), for example, the ground state with the pure Y3 quadrupolar interaction becomes the 
staggered order composed of a sublattice with the state |0) and the other sublattice with the 
state |— ). In this case we expect a large staggered moment induced by the magnetic field. 
The observed staggered field does not support this possibility. 3 Hence in actual PrFe4Pi2, the 
degeneracy of the AFQ should have been lifted already without magnetic field. 

Experimentally, the distortion of ligands observed by X-ray or polarized neutron diffraction 
is [S, 5, <5']-type. 18 ' 19 The dominance of this kind of quadrupolar moment is also indicated by 
the results of polarized neutron scattering, 20 ' 21 where the induced AFM moments is about 2 
times larger for the field direction (0, 0, 1) than for (1, 1, 0). Taking 3z 2 — r 2 -type quadrupolar 
field on one sublattice, the only way to keep the local symmetry at the other sublattice 
is to break the time-reversal symmetry. Namely, if either |+) or |— ) is realized, experiment 
cannot detect the lattice distortion with an O 2 component. However, a large magnetic moment 
should emerge in this case, contrary to the experimental results. On the other hand, if the time- 
reversal symmetry is preserved, the only possible combinations are (|+) + 1— )) and (|+) — | — )). 
In this case, a non-zero x 2 — y 2 -type quadrupolar field should emerge on one of the sublattices. 
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6. Discussion and summary 

In the present work we have proposed a picture for the T3 quadrupole order in PrFe4Pi2 
in terms of the Fi—F^ pseudo-quartet ground state under the symmetry. With the as- 
sumption of vanishingly small energy separation between the singlet and triplet states, the 
low- and the high-field phases can be obtained simultaneously. With the low-lying singlet- 
triplet CEF scheme a level crossing occurs only for the magnetic field direction (111). This is 
the key point for understanding of the appearance of the high-field phase, and the resistivity 
enhancement occurring sharply around this field direction. Introducing ferro-type interactions 
between the dipoles and T@ octupoles, we have reproduced both phase boundaries qualita- 
tively. The magnitude of the octupolar interaction controls whether the high-field phase has 
the same (q = (1,0,0)) or different (q = 0) ordering vector as that in the low-field phase. In 
the latter case we expect macroscopic lattice distortion in the high-field phase. The resultant 
symmetry lower than trigonal should be detected by diffraction measurements. With our set 
of ferro-type interaction parameters the Curie- Weiss temperature ranges from ~ 4.1 to 4.4K, 
which is to be compared with the measured value 0* ~ 3.6K. 

The low-energy CEF level scheme in this compound has not been established, especially 
concerning the separation between the T\ singlet and triplet states. This is related to 
interpretation of the inelastic neutron scattering results that sharp peaks appear only below 
the quadrupolar transition temperature. From our point of view the peak at 1.5meV is ascribed 
to the splitting induced by the AFQ order. Further experiments should identify the origin of 
the inelastic peak, especially by inspection of the temperature dependence near Tq. 

Although the main theme of this paper has been to derive the multiple instabilities of 
the disorderd phase, we have also examined the properties of the quadrupolar order within 
the triplet state by exact two-site and mean-field calculations in the ground state. We have 
found that the AFQ ordering of the T3 quadrupoles in the T/j triplet results in a degeneracy 
with respect to different quadrupolar patterns. Other mechanism beyond the T% quadrupole 
interaction is necessary to resolve the degeneracy of the quadrupoles in the triplet. At present, 
we cannot exclude the possibility that a tiny x 2 — y 2 -type distortion is actually realized. 

In order to discuss the lattice distortion microscopically, we have to consider the coupling 
of the T3 quadrupoles to the lattice of Fe and/or P ions. Due to the bcc lattice structure, 
the Fe ion lattice is shared between Pr ions located on different sublattices. Because of the 
" non-collinear" antiferro-quadrupolar order, the distortion of the Fe ion lattice required by 
the quadrupolar moment on one site is not the optimal one for the other site. 
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